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ABSTRACT 

We compute the linear power spectrum of cosmic string sedeed fluctuations in 
the context of neutrinos with a strong self-interaction and show that it is very similar 
to that obtained in the context of 'normal' neutrinos. We compare our results with 
observational data and show that for any value of the cosmological parameters h and 
Qq the interacting hot dark matter power spectrum requires a large scale dependent 
biasing parameter. 
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1 INTRODUCTION 

Observations on all scales of cosmological interest show that 
most of the matter in the universe is in some form of un- 
seen matter. On the other hand, nucleosynthesis constraints 
seem to imply that most of this dark matter is non baryonic. 
Gravitinos, axions and neutrinos for example, are only a few 
among the many candidates to constitute the dark matter. 

In the standard scenario it is usually assumed that 
during the important epochs for structure formation, the 
dark matter particles candidates interact only gravitation- 
ally. However one could assume that this is not the case. 
In particular it is possible that light neutrino-like particles 
could have some other kind of self-interaction during the low 
energy epochs when structure formation takes place (see, for 
instance, Raffelt and Silk 1987; Carlson, Machacek and Hall 
1992; Gradwohl and Frieman 1992; Machacek 1994; Laix, 
Scherrer and Schaefer 1995). 

In recent work Atrio-Barandela and Davidson (1997), 
have considered a light (~ 30 eV) self- interacting neutrino- 
like particle and discussed the possibility that the dark mat- 
ter in the universe could be constituted by this kind of parti- 
cles. They determined the linear power spectrum of density 
fluctuations generated by the present time in the context of 
primordial gaussian fluctuations and concluded that galaxy 
sized density perturbations could survive. 

In what concerns the estimates for the bounds on the 
neutrino-neutrino cross section very little is known. Within 
the framework of the Standard Model of Particle Physics, 
no experimental bounds have been set, and the best known 
lower limit on the mean free path of the neutrinos is based on 
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observations of the supernova SN 1987 A. In fact, the neutri- 
nos produced by SN 1987^4 in the Large Magellanic Cloud, 
reached the earth after crossing the relic neutrinos of the 
cosmic background. Since the number of detected neutrinos 
agrees with the theoretical expected number from supernova 
models, the present mean free path of the neutrinos Ao, is 
not shorter than the distance to SN 1987 A, which is of the 
order of 2 x 10 23 cm; for a complete analysis, see Kolb and 
Turner, 1987; (see also Nussinov and Roncadelli, 1983). In 
this article we will assume that the co-moving value of A is 
much smaller than any scale of cosmological interest dur- 
ing the time where most of the structures we observe in the 
universe today are generated. In this case one can treat the 
neutrino component with the hydrodynamic or fluid approx- 
imation. 

Here we extend the work of Barandela and Davidson 
(1997) to the case of non-gaussian density perturbations in- 
duced by cosmic strings. Cosmic strings, along with other 
topological defects, may be produced as a result of symme- 
try breaking phase transitions in the early universe (Kibble 
1976). Although some topological defects are ruled out by 
observations, others like cosmic strings could be the seeds 
which generated the large scale structures we observe in the 
universe today. They also have other interesting cosmolog- 
ical consequences (see for example Vilenkin 1985; Colombi 
1993; Vilenkin and Shellard 1994; Avelino 1996). Among 
these there is a non-gaussian signature in the microwave 
background (on small scales), the production of double im- 
ages and a primordial background of gravitational waves or 
other radiation. 

The outline of this paper is as follows. In section 2 we 
start by describing the formalism we employ in order to 
study the growth of density fluctuations in the context of 
neutrinos with a strong self-interaction. We introduce com- 
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pensation in section 3 and in section 4 we describe the semi- 
analitic model we employ in order to compute the power 
spectrum of cosmic string seeded density fluctuations. In 
section 5 we describe and discuss the results. We conclude 
in section 6. 

In this article employ fundamental units with % = c = 
fcs = 1. We consider a self-interacting neutrino with mass 
m„ = 93Qoh 2 such that fio = f2„. 



2 SUBSEQUENT PERTURBATIONS 

In a flat FRW universe with no cosmological constant and 
containing both CDM and radiation fluids, the scale factor 
a may be written in terms of the conformal time 77 as 



a(rj) = a e 



2(V5-l)-2- + (3- 

Veq 



(1) 



where a eq and r\ eq are respectively the expansion factor and 
conformal time at matter-radiation equality. If the dark mat- 
ter is hot but the transition from the relativistic to the non- 
relativistic regime occurs deep in the radiation era eqn. (|l]) 
is still a good approximation. 

The evolution of radiation and non-relativistic interact- 
ing hot dark matter (IHDM) density fluctuations in the syn- 
chronous gauge is given by 
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where Q a p is the energy-momentum tensor of the external 
source, c s is the adiabatic sound speed of the interacting hot 
dark matter and ilh and il r express the densities in radiation 
and IHDM as fractions of the critical density. 

We are implicitly assuming that the strength of the 
neutrino-neutrino coupling is large enough for it to be a good 
approximation to treat the neutrino component as a perfect 
fluid. We also assume the radiation behaves as a fluid; this is 
certainly valid before recombination for the length scales of 
interest, well above the radiation damping scale. Although 
this is not true after recombination, for Qoh 2 not too small 
the radiation no longer dominates the dynamics after last 
scattering. Consequently, we expect that even then equation 
(^) remains a valid approximation. 

The solution to the system of equations with initial 
conditions Sh — 8 r — Sh — S r = can be written in terms of 
Green functions as 

5f, r (x^)=47r r dr,' [ d 3 x'g h - r (X; V ,r,')e + (x', V '), (4) 
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Here, O+ = 0oo + &u and X = x — x'|. The upper index 
'S' indicates that these are the 'subsequent' fluctuations, ac- 
cording to the notation of Veeraraghavan & Stebbins (1990), 
to be distinguished from the 'initial' fluctuations. In Fourier 
space, the Green function Q ' r obey the same equations as 

&h,r- 



Figure 1. Plot of T, fi p (in units of m), and c s as a function of 
the scale factor a normalized to unity at r; eq . 



gh + V _ I (1) {nh g* + 2 n r g r ) + c 2 s k 2 g h 

a 2 \a/ 
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The initial conditions at r\i are g h = 1, <J r = 4/3 and g h = 
g T = 0. Hence, given the evolution of the source stress- 
energy, it is possible to numerically compute the resulting 
density fluctuations. 

We are mostly interested in the interacting IHDM in- 
homogeneities at late times in the matter era. In the limit 
T) S> rj eq , the Green functions are dominated by the growing 
mode, oc a/a eq . Hence, the function we would like to solve 
for is 



T(k;r]i 



lim ^lg*(k,n,r)i). 



(8) 



This function will be used later to construct the transfer 
function for the power spectrum of IHDM density fluctua- 
tions induced by cosmic strings. 



2.1 The IHDM sound speed 

The adiabatic sound speed of the self-interacting neutrinos 
c s , is given by 



(9) 



where p stands for pressure, p is the mass density of the 
neutrinos and the subscript S on the right hand side means 
that the entropy per particle S, is constant. 

The neutrinos have Fermi-Dirac distribution function. 
Therefore, the neutrino number density n, mass density p 
and pressure p are given respectively by: 
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E (E 
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dE, (10) 
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fp 2 + to 2 ) 1 2 is the energy, p is the momentum, 



where E — 

m is the mass, T is the temperature, p p = p P (T) is the 
chemical potential, A = g/(2n) 2 , and g is the spin degener- 
acy factor of the neutrinos. Using the constraint 

Mp ) = l5(3)' (13) 

where S is the entropy per neutrino, we numerically calcu- 
lated n, p, p and functions of T. We then obtained the 
adiabatic sound speed c 2 = (dp/dp)s, and fit the obtained 
result by: 



S = 



I ( P + > 
TV n 



Cs(a.) = — 
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-1/2 
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where a = 0.19 is the best fit parameter. In Fig. 1 we rep- 
resent the numerical plots of T and p p (in units of m), and 
c s , as functions of the scale factor a, normalized to unity at 
the epoch of equality between matter and radiation. 



Figure 2. Comparison between the CMB-normalized linear 
power spectrum of cosmic string seeded fluactuations with CDM 
(dot-dashed line), IHDM (solid line) and HDM (doted line) for 
Q.q = 1 and h = 1. 



3 COMPENSATION 

The linear perturbations induced by cosmic strings, are the 
sum of initial and subsequent perturbations: 

Sh{k;ri) = 5i(k;rj) + 6h(k;rj) 



= 4n(l + z eq ) / d n 'T{k-ri')Q + {k- n '). (15) 

The transfer function for the subsequent perturbations, 
those generated actively, was obtained in the previous sec- 
tion. To include compensation for the initial perturbations, 
8^, we make the substitution: 



T(k; V )^ (l + (k c /k) 2 ) T(k; V ) 



(16) 



where k c is a long-wavelength cut-off at the compensation 
scale. This results from the fact that local physical processes 
cannot produce perturbations on scales much larger than 
the horizon size. Consequently, the power spectrum of den- 
sity perturbations is bounded by a k 4 spectrum on large 
scales (k << k c ), assuming that background fluctuations 
are uncorrelated for points which have never been in causal 
contact. In the I-model of strings, which most resembles 
the Bennett-Bouchet (1990) and Allen-Shellard simulations 
(1990), k c = 271-/77. 



4 POWER SPECTRUM 

The analytic expression of Albrecht and Stebbins (1992) for 
the power spectrum of density perturbations induced by cos- 
mic strings in a fl = 1 FRW universe with no cosmological 
constant is 

P(k) = 16TT 2 (l + z eq ) 2 (Gp) 2 



dr,'F[kZ(v')/a(v'W(k, V ')f 



(17) 



where the "p + 3p" part of the string stress energy tensor, 
is modeled by T given by 



FU/a] = J^/3 2 E|!(l + 2(fc X /a 



(18) 



In these equations, ryo is the conformal time today, rji is the 
conformal time when the string network was formed, \ is 
the typical curvature scale of the wakes, 



i = { Pao /ny, f3 = (v 2 y, (19) 

+ (20) 

where is the energy density in long strings, p, is the string 
mass per unit length, v is the microscopic velocity of the 
string, Pb is the macroscopic bulk velocity of the string, 75 = 
(1 — /9 2 ) 5 , and p r is the renormalized mass-per-unit-length. 
We note that E/3 2 ~ 1 and x/ a ~ 2 U a ~ both in the 
matter and radiation eras (Bennett and Bouchet 1990; Allen 
and Shellard 1990) and consequently the structure function 
T is always well approximated by 

1 



s 



1 + ± (M 2 



(21) 



We substitute this in equation (^) in order to calculate the 
IHDM power spectrum. 



5 RESULTS 

In Fig. 2 we compare the power spectrum of the perturba- 
tions generated by cosmic strings for both CDM (dot-dashed 
line), HDM (doted line) and IHDM (solid line) for ft = 1 and 
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lems as the standard HDM scenario for cosmic strings. We 
shall discuss these problems next when we investigate other 
choices for the cosmological parameters. 



5.1 Open models 

The generalization of these results for open models of struc- 
ture formation can be made using the same technics em- 
ployed by Avelino, Martins and Caldwell (1997). In Fig. 3 
we compare the power spectrum of the perturbations gen- 
erated by cosmic strings in the context of IHDM with the 
Peacock and Dodds (1994) linear power spectrum reconstru- 
cion inferred from various galaxy surveys. 

The IHDM power spectrum is again normalized to the 
COBE-DMR observations (Allen el al., 1997). We see that 
in an open universe it is possible to obtain a better fit to the 
shape of the Peacock and Dodds linear power spectrum on 
large scales. However, in this case perturbations are damped 
on too large scales which is in conflict with the observational 
data. In a flat universe with a cosmological constant the 
power spectrum requires a slightly lower biasing than for 
an open universe with the same matter density (Avelino, 
Caldwell and Martins, 1997). However, because the shape of 
the power spectrum density fluctuations induced by cosmic 
strings remains the same this model is also in conflict with 
observations. 



Figure 3. Comparison between the CMB-normalizcd linear 
power spectrum of cosmic string seeded fluactuations with IHDM 
(solid line) for h = 0.7 and Cl = 1.0,0.4,0.2. For H < 1 the 
reconstructed linear power spectrum has been rescaled by f2~°' 3 
(Peacock and Dodds, 1994). 



6 CONCLUSIONS 

In this paper we calculated the power spectrum of density 
fluctuations induced by cosmic strings in the contex of neu- 
trinos with strong self-interactions. We concluded that be- 
cause gravitational instability can only be established on 
scales smaller than the Jeans scale, small scale power is re- 
moved relative to the CDM case. In opposition to the case 
studied by Barandela and Davidson (1997), we note the ab- 
sence of oscillations, which is due to the incoherent nature 
of the cosmic strings source. In fact, the results obtained 
for HDM and IHDM are very similar, the IHDM spectrum 
having a slightly larger amplitude on small scales. The gen- 
eralization of these results shows that in an open universe 
or a flat universe with a cosmological constant it is possi- 
ble to decrease the biasing on large scales by a small factor. 
However, in this case perturbations are damped on too large 
scales. Hence, we conclude that it seems difficult to reconcile 
these results with observations unless there is some physi- 
cal mechanism which is capable of generating a large, scale 
dependent biasing. 



Qo = 1. The curves are normalized to the COBE-DMR ob- 
servations (Allen el al., 1997). We can see that the HDM 
and IHDM power spectra are very similar, the IHDM power 
spectrum having a slightly larger amplitude on small scales. 
The absence of oscillations reflects the incoherent nature of 
string perturbations which is implicit in the way we add the 
perturbations generated by the cosmic strings at different 
times to calculate the power spectrum. Hence, it seems that 
the IHDM cosmic string scenario suffers from similar prob- 
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